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The LDA+DMFT (local density approximation combined with dynamical mean-field theory) com- 
putation scheme has been used to study spectral and magnetic properties of FeSi and Fei-xCo^Si. 
Having compared different models we conclude that a correlated band insulator scenario in contrast 
to Kondo insulator model agrees with FeSi band structure as well as experimental data. Coulomb 
correlation effects lead to band narrowing of the states near the Fermi level with mass renormaliza- 
tion parameter m* 2 in agreement with the results of angle-resolved photoemission spectroscopy 
(ARPES). Temperature dependence of spectral functions and magnetic susceptibility calculated in 
DMFT reproduces transition from nonmagnetic semiconductor to metal with local magnetic mo- 
ments observed experimentally Cobalt doping leads to ferromagnetism that has itinerant nature 
and can be successfully described by LDA+DMFT method. 

PACS numbers: 71.27.+a, 71.10.-w 



I. INTRODUCTION 

The narrow-gap semiconductor FeSi demonstrates an 
interesting interplay between magnetic and electronic 
properties. Magnetic susceptibility temperature depen- 
dence shows maximum at 500 K and Curie- Wiess behav- 
ior at higher temperatures.^ While intrinsic magnetic sus- 
ceptibility vanishes below 50 K, FeSi does not show any 
sign of spin ordering down to the lowest temperatures.^ 
Photoemission^ and optical experiments show an energy 
gap of about 60 meV at low temperatures that is gradu- 
ally filled with temperature increase.^ A resistivity tem- 
perature dependence shows transition from a narrow gap 
semiconductor to a bad metal. ^ Cobalt doping results in 
ferromagnetic metal state for Fei_a;Coa^Si.^ 

Several models have been suggested to explain the un- 
usual temperature dependence of the FeSi physical prop- 
erties ranging from spin fluctuations^i to phenomenolog- 
ical models assuming two narrow (i-bands in the vicinity 
of the band gap.^^^ Such density of states (DOS) mod- 
els are similar to a Kondo insulator description, and due 
to the striking similarities in the physical properties of 
these Kondo insulators it was claimed that FeSi is the 
first Kondo insulator containing no /-electrons ^iSiiii 

Historically, the first FeSi model was proposed by Jac- 
carino et al.^ To describe the unusual magnetic suscep- 
tibility and specific heat a model DOS was proposed 
with extreme narrow band peaks around a small en- 
ergy gap. Despite the fact that the model fitting results 
were in good agreement with experimental data for the 
susceptibility this picture contradicts a band structure 
calculation^^ where no unphysical narrow bands were ob- 
tained. 

The next model proposed to describe electronic and 



magnetic properties of FeSi was the Kondo insulator 
model. The motivation of applying the Kondo model 
to explain physical properties of FeSi was that the spin- 
fluctuation spectra of CeNiSn and FeSi are similar. Ac- 
cording to this model, a set of localized atomic-like elec- 
tron levels interact with a wide itinerant band. The in- 
sulating state scaled by a Kondo temperature T^ is the 
result of a weak hybridization between the localized and 
itinerant bands. The implementation of the Kondo in- 
sulator model to FeSi compound is however questionable 
since the band structure of FeSi shows a strong hybridiza- 
tion between the Fe-3(i and Si-p states. Also the mag- 
netic interactions in FeSi are essentially ferromagnetic 
and not antiferromagnetic, as it would be expected from 
an RKKY-picture. For all these reasons it is desirable 
to have a microscopic model based on the realistic band 
structure which can reproduce experimental data. 

First-principle band structure calculations of Mattheis 
and Hammann^^ have shown that energy gap value of 0.1 
eV is essentially smaller than the width of the band above 
the Fermi level (~ 0.5 eV). Band structures analysis led 
the authors to conclusion that hybridization between Fe- 
d and Si-p states is very strong. These results do not 
support the Jaccarino's model as well as the Kondo sce- 
nario. 

Recently, another model of a correlated band insula- 
tor (CBIM) was proposed by Kunes and Anisimov^i^ Us- 
ing the Dynamical Mean-Field Theory the authors have 
taken into account local dynamical correlations for a 
gapped LDA spectral function of FeSb2 that experimen- 
tally demonstrates transition with temperature increase 
from nonmagnetic semiconductor to metal with local mo- 
ments similar to FeSi. It was found that the DMFT en- 
ergy gap is reduced due to correlation effects by a factor 
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of two from its LDA value. Within the CBIM picture 
bands above and below the energy gap are formed by 
nonlocal bonding and antibonding orbital combinations. 
The on-site Coulomb interaction leads to a competition 
between localization and formation of nonlocal bonds. 
The DMFT calculations allowed to reproduce success- 
fully temperature dependence of magnetic susceptibility, 
resistivity and optical conductivity experimentally ob- 
served in FeSb2. 

Recent angle-resolved photoemission experiments^^ 
have proved validity of the correlated band insulator sce- 
nario for FeSi. The authors have found an effective band- 
mass renormalization m*/m?^2. As a consequence the 
LDA band gap of 100 meV is renormalized by a factor 
of 0.5 to about 60 meV in agreement with other experi- 
mental data. 

We summarized the parameter values of the models 
used to describe FeSi in Table [H The basic difference 
between Kondo insulator and correlated band insulator 
(CBIM) models is a ratio between energy gap value Egap 
and a width W of the bands around the gap. While 
Kondo insulator model requires W/Egap <<1 to give 
good results, CBIM uses W/Egap >>1 in agreement with 
FeSi realistic band structure. 



TABLE I: Classification of FeSi models. 



Jaccarino's modeii W « Egap (W ^ K, Egap = 1520 K) 
Kondo 

insulator model^ W = Egap/2 (W = 500 K, Egap = 1000 K) 
CBIM (this work) W » Egap (W = 5000 K, Egap = 1000 K ) 



II. FeSi 

A. Band effects 

The crystal structure of FeSi corresponds to four for- 
mula units in unit cell and its band structure is rather 
complicated for analysis using a simple model. To pro- 
vide a better understanding of FeSi band structure near 
the Fermi level Mattheiss and Hamann^^ have proposed 
to consider the closely related phase with the symmetry 
of the rocksalt structure. The latter phase can be ob- 
tained from original crystal structure by relatively small 
shift of atomic positions. They have shown that the ori- 
gin of the FeSi energy gap can be traced to a pseudogap 
that is present in the rocksalt phase. This result is a 
starting point of our investigation that is aimed to con- 
struct a minimal realistic band insulator model which 
captures main features of LDA electronic spectrum of 
FeSi. 




FIG. 1: TB-LMTO energy band of FeSi obtained with dif- 
ferent sets of atomic position parameters u(Fe) and u(Si). 
(a) Structure with u(Fe)=0.25 and u(Si)=0.75 corresponds to 
nonprimitive rocksalt structure that contains four FeSi for- 
mula units. Figures (b) and (c) are energy band results for 
two transitional structures with u(Fe)=0.23,u(Si)=0.76 and 
u(Fe)=0.19,u(Si)=0.79, respectively. Figure (d) corresponds 
to the real FeSi structure. Symmetry lines are chosen accord- 
ing to Ref. il2| 

As first step we have studied effects of atomic posi- 
tions shift on the band picture of FeSi. For that pur- 
pose we have carried out calculations using the Tight 
Binding Linear-Muffin-Tin- Orbital Atomic Sphere Ap- 
proximation (TB-LMTO-ASA) method with conven- 



This paper is aimed to provide microscopic analysis of 
hybridization and correlation processes in FeSi. First we 
propose a simple band insulator model that captures the 
essential hybridization effects between Fe-3d and Si-3p 
states in FeSi (Section III Ap . This model demonstrates 
an energy gap that is five times narrower than the width 
of bands. In order to investigate correlation effects we de- 
fine a density of states model based on LDA calculations 
for FeSi and solve it using the Dynamical Mean-Field 
Theory (Section III Bp . A strong renormalization of DOS 
near the Fermi level was found in the DMFT calculations 
that is in good agreement with a band narrowing ob- 
served in the recent ARPES experiments. Temperature 
increase in DMFT calculations results in transition from 
nonmagnetic insulator to a bad metal with local moments 
in agreement with experimental data. Having encour- 
aged by these results we have investigated a doped FeSi 
model to study the magnetic properties of Fei-a^Coa^Si 
alloys (Section We have demonstrated that a good 
agreement between calculated and experimental values 
for magnetization M as a function of doping x can be 
achieved using the LDA+DMFT method. 
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tional local- density approximation (LDA)^^ for real (sim- 
ple cubic) structure (atomic positions u(Fe) =0.1358 
and u(Si) =0.844), rocksalt (face-centered-cubic) phase 
(u(Fe)=0.25 and u(Si)=0.75) and for two model struc- 
tures with atomic positions intermediate between those 
values (Fig. [T]). The obtained bands agree with LAPW 
results for real and rocksalt phases. The fee band struc- 
ture contains several electrons and holes pockets near G 
and X points in the Brillouin zone that corresponds to 
a pseudogap in the energy spectrum in contrast to the 
real gap in the simple cubic structure. For intermediate 
structure with atomic position parameters u(Fe)=0.19 
and u(Si)=0.79, the energy gap is open along all symme- 
try lines of the Brillouin zone with the exception of GM . 
One can see that there is a complete energy gap in the 
real FeSi structure. This result agrees with a conclusion 
of Mattheiss and Hamann^^ that the energy gap state 
of FeSi results from the distortion of the pseudogapped 
rocksalt structure. 

It is naturally to expect that the features of the FeSi 
electronic spectrum (such as the narrow gap and the peak 
above the Fermi level) are provided by Ye-Zd and Si-3s, 
3p states J^ii^ To investigate the origin of the pseudogap 
state in the fee phase we artificially scaled the hybridiza- 
tion strength between Si-3s, 3p and Fe-3(i states. To do 
so off-diagonal Hamiltonian elements which describe the 
hybridization between Fe and Si atoms were multiplied 
by a coefficient < a < 1. These results are presented 
in Fig. [2 

To simplify our analysis of the band structures one 
can consider dispersion curves along LG symmetry line 
where complete energy gap opens due to hybridization 
between Ye-Zd and Si-35, 3p states. Without hybridiza- 
tion (a = 0) the wide silicon band crosses three narrow 
3(i-bands of iron for /c-vectors along LG direction. At 
increasing hybridization strength (Fig. [2)3 and[2]3) two of 
the iron bands strongly interact with silicon states while 
one of them remains practically unaffected hy d — p hy- 
bridization. The resulting dispersion curves of the fully 
hybridized system are presented in Fig. [2] (d). It is clear 
that the pseudogapped state (Fig. [2]i) originates from 
a strong hybridization of the wide silicon band (~ 10 
eV) with the relatively narrow iron band (~ 1 eV). The 
narrow gap opens between antibonding hybridized band 
above and unhybridized band below the Fermi level. 

While qualitatively this situation resembles the Kondo 
insulator picture with weakly hybridized wide and narrow 
bands, there are two essential quantitative differences. 
First of all the "narrow" band having width of an order 
of magnitude smaller than the "wide" band is still too 
wide in absolute value of about 1 eV. At second d — p 
hybridization is so strong that "pure hybridization" gap 
in Fig. [2]i is larger than 2 eV. The small value of the gap 
is not due io d — p hybridization weakness but happens 
between hybridized and unhybridized bands. 

Another argument against the Kondo scenario can be 
found from analysis of the two bands (Fig. [T]i) form- 
ing a well-separated narrow peak above the Fermi level 
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FIG. 2: Band structures for the different degree of the hy- 
bridization strength between Fe-3d and Si- 3s, 3p states in fee 
crystal structure with one formula unit in unit cell. Figures 
(a), (b), (c) and (d) are the energy spectra for hybridization 
scaling parameter a=0; 0.1; 0.5 and 1, respectively. Symme- 
try lines are chosen according to Ref. 



(Fig. [6]). In Kondo system these bands should be asso- 
ciated with strongly localized atomic orbitals of an iron 
atom. To check this we have calculated Wannier func- 
tions for two bands above the Fermi level using a projec- 
tion procedure^ii A spatial distribution for one of the cal- 
culated Wannier functions is shown in Fig. [3l The result- 
ing Wannier functions correspond to a complex combina- 
tion of 3(i-states of iron and 3s-, 3p-states of silicon and 
are spread over whole unit cell of FeSi containing eight 
atoms. This picture is very far from localized atomic 
orbital needed for the Kondo scenario and supports the 
band insulator model. 

We built an effective microscopic model (Fig. |4]) for LG 
direction, that contains a minimal set of orbitals and re- 
produces the gapped state. The model silicon and model 
iron are described with one and two orbitals, respectively. 
The model Hamiltonian is given by 



H 
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where ti, ^2, ts and ^4 are hoppings between model iron 
and silicon orbitals presented in FigHl We have esti- 
mated the hopping parameters of the LG model using 




FIG. 3: (Color online) Wannier function corresponding to 
a narrow band above the Fermi level (see Fig. [TJi). The 
Wannier function is centered at a 3d orbital of iron atom 
(red spheres). We found about 40% of the electron density 
at the central atom and its iron neighbours. Green spheres 
correspond to silicon atoms. 



the real band structure of the full Hamiltonian presented 
in Fig. [21 The obtained hopping integrals are ti=2.0 
eV, t2=3.0 eV, t3=-1.2 eV, t4=-0.25 eV. These values are 
much larger than those calculated for localized systems^ 
and are far outside of the values range needed for the 
Kondo insulator scenario. 

The calculated model band structures with {a=l) and 
without {a=0) hybridization are presented in Fig. [5l One 
can see that our model results are in good agreement with 
those obtained from calculations with full Hamiltonian 
(Fig. [2]). The substantial Fe-Si hybridization leads to 
splitting into a lower bonding and an upper antibonding 
bands with a non-bonding band in between. We have 
obtained a small gap semiconducting ground state with 
the energy gap value of 0.8 eV that is much smaller than 
the narrowest band width of 4.8 eV near the Fermi level. 
One should note that the small value of the energy gap 
is not a consequence of a weak Fe-Si hybridization and 
the minimal realistic band structure model of FeSi differs 




FIG. 4: (Color online) One-dimensional effective model. 
Large and small circles correspond to iron and silicon atoms 
in the fee phase of FeSi. 
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FIG. 5: Model band structures. Lower and upper figures 
correspond to the models with and without hybridization, re- 
spectively. 



qualitatively from the Kondo insulator regime. 

The proposed band model can be used to study 
Coulomb correlation problem within statical or dynami- 
cal mean- field approaches. However, as we will show be- 
low, it is possible to define an effective density of states 
model derived from the full DOS obtained in LDA cal- 
culation for the real crystal structure. Solution of this 
model by DMFT allows us to describe anomalous physi- 
cal properties of FeSi. 



B. Correlation effects 

In this section we study the effects of the Coulomb cor- 
relation on the electronic structure and physical proper- 
ties of FeSi using the Dynamical Mean-Field Theory. As 
an input DMFT requires a non-interacted Hamiltonian 
or a density of states (DOS).^^ The essential features of 
FeSi density of states obtained in LDA calculations (Fig. 
[6]) are small energy gap ^100 meV and narrow (^^0.5 eV) 
peak above the Fermi level containing 0.5 electrons per 
spin per Fe atom. We have defined the model density of 
states (see filled area in Fig. [6]) by cutting from the entire 
DOS area around Fermi energy with an integral equal to 
1 electron per spin per Fe atom. The model density of 
states contains the main features of the FeSi spectrum, 
namely the energy gap (0.1 eV) and the narrow peak 
above the Fermi level (0.5 eV). 

As impurity solvers of the DMFT problem we have 
used a Quantum Monte-Carlo method with Hirsch-Fye 
algorithm^^ (QMC-HF), a Continuous-time Quantum 
Monte- Carlo method with interaction expansion^ ^ (CT- 
QMC) and an exact diagonalization temperature depen- 
dent Lanczos^^ approach. 

The value of the on-site Coulomb interaction param- 
eter U was chosen to be 1 eV that is close to the value 
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used in Ref. 14. This rather smah value can be justified 
by an effective screening of d—d Coulomb interaction due 
to strong Fe-ligand hybridization as it was demonstrated 
in constrain DFT calculation of the Coulomb interaction 
parameter U for LaOFeAs.^^ 
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FIG. 6: (Color online) Model (filled area) and full (line) den- 
sities of states. 

The resulting paramagnetic densities of states at T = 
232 K calculated in DMFT using various impurity solvers 
are presented in Fig. H One can see that all methods re- 
sult in a pseudogapped state with strongly renormalized 
density of states near the Fermi level. There are satel- 
lites at ±17/ 2 which correspond to the lower and upper 
Hubbard bands. 

The energy area around the Fermi level corresponds 
to quasiparticle states that are usually described as non- 
interacting bands renormalized by Coulomb correlations. 
The renormalization process can be understood by a low- 
frequency analysis of the dynamical mean- field equations. 
The lattice Green function is given by 



Gicu) = J2icu - Eicu) - e(k))- 



(1) 
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FIG. 7: (Color online) Spectral functions obtained from 
DMFT calculations using QMC-HF (blue dashed line), CT- 
QMC (black solid line) and exact diagonalization (red solid 
line) techniques at T = 232 K. Orange filled area corresponds 
to LDA density of states. 
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FIG. 8: (Color online) Density of states obtained from QMC 
calculations at T= 386 K (black line) and 96 K (orange line). 



where S is a self-energy and e(k) is a LDA spectrum. We 
expand the real part of the self-energy in the vicinity of 
the Fermi energy leaving only linear term: 

Re^ico) « i?eS(0) + oj^^^^l^o- (2) 
Then Green function for specific wave vector k is 

GkM«(u;m*-e(k))-i, (3) 

where m* is the effective band mass renormalization pa- 
rameter 



m* = l v^U=o. 

duj 

Equation ([3|) can be rewritten as 



w - e(k) ' 



(4) 



(5) 



where Z = 1/m* is a quasiparticle weight and e(k) = 
e(k)/m* is a new band dispersion. Therefore, the 
Coulomb correlation renormalization for quasiparticle 
states near the Fermi level results in a band structure 
narrowing by a factor of m* and a reduction of the corre- 
sponding spectral weight by a factor of Z with the rest of 
the spectral weight transfered to upper and lower Hub- 
bard bands at ±/7/2. 

In our DMFT calculations the effective band mass 
renormalization parameter m* ~ 2 was found in good 
agreement with the results of angle-resolved photoemis- 
sion spectroscopy (ARPES)^^ showing band narrowing 
by a factor of two in comparison with bands calculated 
in LDA. 

Experimentally FeSi demonstrates transition with 
temperature increase from narrow gap semiconductor to 
bad metal*^ We have performed DMFT calculations at 
different temperatures and such an experimentally ob- 
served transition was successfully reproduced theoreti- 
cally. The DMFT spectral functions calculated for differ- 
ent temperatures are presented in Fig. [HI The energy gap 
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FIG. 9: (Color online) The convolution for FeSi model at 386 
K (black line), 232 K (red line), 116 K (brown line) and 96 
K (orange line). Inset: The experimentally observed optical 
conductivity taken from Ref. 



of about 50 meV found at T= 96 K agrees with experi- 
mental results of resistivity measurements which indicate 
a charge gap of about 60 meV.^"^ At T= 386 K the gap 
area is nearly completely filled by spectral weight trans- 
fer from sharp peaks near the gap resulting in the spec- 
tral function corresponding to a bad metal. The energy 
gap value obtained in the DMFT calculation is two times 
smaller than corresponding value from LDA band struc- 
ture calculations (0.1 eV). That agrees very well with the 
effective band mass renormalization parameter m* ~ 2 
Eq.dlj) obtained in our DMFT calculations. 

Optical spectroscopy experiments show gradual fill- 
ing of the low-temperature energy gap with tempera- 
ture increase till complete gap disappearance above room 
temperature^^ We have estimated the optical conductiv- 
ity by a spectral function convolution using the following 
expression 

a{u,T) = - /de7V(e)7V(e + c^)[l-/(e + c^,T)], (6) 

where N{e) is a spectral function obtained from DMFT 
calculations (Fig. [7j) and f{e^T) is the Fermi distri- 
bution function. The calculated convolution together 
with experimental data for temperature dependent op- 
tical conductivity^ are presented in Fig. [9l 

At low temperatures a well pronounced energy gap of 
about 0.05 eV can be observed in both experimental and 
theoretical curves. With temperature increase optical 
conductivity increases at the energies below 0.05 eV and 
finally at T= 386 K there is no any trace of the gap in 
the theoretical curve in good agreement with experimen- 
tal data. 

FeSi displays an unusual crossover in the vicinity of the 
room temperature from a singlet semiconducting ground 
state with a narrow band gap to a metal with an en- 
hanced spin susceptibility and Curie- Weiss temperature 
dependence.^ All previous attempts to explain this be- 
havior were based on the models assuming extremely nar- 
row (< lOOOK) peaks at the energy gap edges in DOS, 



while LDA calculations gave band width nearly an order 
of magnitude larger than that value. Our DMFT cal- 
culations demonstrate that the correlated band insulator 
model with realistic DOS can reproduce anomalous tem- 
perature dependence of the magnetic susceptibility for 
FeSi. We have computed x(^) a ratio 



X(T) = 



M 



(7) 



where his a small uniform external magnetic field and M 
is an induced magnetization of the system. The compar- 
ison of the experimental and calculated magnetic suscep- 
tibilities is presented in Fig. [TOl The calculated spin sus- 
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FIG. 10: (Color online) Spin susceptibility x{T) ( in //|/ eV) 
from LDA+DMFT calculations (orange solid line) and from 
experiment^ (blue dashed line). 

ceptibility increases exponentially from T=0 and reaches 
maximum at 600 K what is in good agreement with the 
experimental temperature of 500 K. However, the cal- 
culated absolute values of the magnetic susceptibility at 
T > 300 K disagrees with experimentally observed x(^)- 
Such a disagreement was observed in the previous theo- 
retical investigation.^^ Varying the Coulomb interaction 
U the authors of Ref. 25 were able to fit a theoretical 
curve to the experimental one. 

In the present investigation the same model DOS (Fig. 
[6]) and Coulomb interaction parameter U value was used 
to reproduce successfully band narrowing observed in 
ARPES/^ optical conductivity^ and temperature depen- 
dence of magnetic susceptibility^ 



III. Fei_xCoxSi 

The correct theoretical description of magnetic proper- 
ties for transition metal monosilicides such as for exam- 
ple MnSi and Fei-a^Co^^Si alloys presents a longstanding 
problem. There has been a considerable amount of exper- 
imental and theoretical work on MnSi and Fei-a^Co^^Si, 
regarding their structural, magnetic and electronic prop- 
erties. However, at the moment there is no satisfac- 
tory first-principles description of monosilicides magnetic 
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properties. For instance, in case of MnSi it was found 
that the experimental value of the magnetic moment is 
about 0.4 /JjB' Different first-principles calculations based 
on the density functional theory give much larger mag- 
netic moment value of 1 jub^^ 

Fei-a^Coa^Si alloys are magnetic for almost all of the 
intermediate concentration regimes, ^^^^^^^ while the end 
compounds FeSi and CoSi are nonmagnetic, the latter 
being a diamagnetic semimetal. Fei-a^Co^^Si system is 
also interesting for scientists due to the promising prop- 
erties for spintronic device applications. For instance, 
in paper— the authors have reported the discovery of a 
large anomalous Hall effect for Fei-a^Co^^Si. They have 
demonstrated that the large effect is most likely intrinsic 
- derived from the band-structure effects rather than due 
to impurity scattering. They have proposed to consider 
the transition metal monosilicides as potential alterna- 
tives to the (GaMn)As and (GaMn)N which are the most 
popular materials for spintronics. 

From theoretical side no calculation reported so far 
seems to reproduce correctly both magnetic moment 
value and Curie temperature of Fei-a^Co^^Si system. To 
model Fei_a^Coa^Si alloys the authors of paper^^ have 
used the full potential linearized augmented plane-wave 
method in combination with a virtual crystal approxi- 
mation (VGA) as well as with a supercell approach. The 
resulting magnetic moments agree with experimental val- 
ues only for x < 0.25. Having supposed an important 
role of ordering and segregation they have simulated sev- 
eral alloy configurations for concentrations x < 0.5. The 
weighted average magnetic moment which was calculated 
through Boltzmann distribution is in reasonable agree- 
ment with experimental values. However, there are no 
experimental results which support the segregation and 
ordering phenomena. Moreover, the results of neutron 
measurements^^ have demonstrated a random distribu- 
tion of transition metals in Feo.5Coo.5Si. 

To simulate randomly distributed Fei-a^Co^^Si alloys 
the authors of Ref. [H have used the combination of an 
exact muffin-tin orbitals method and a coherent potential 
approximation. They have found an extreme sensitivity 
of magnetic properties to the internal structure parame- 
ters and lattice constant. However, the calculated mag- 
netic moments at concentrations x > 0.3 still disagreed 
with those experimentally observed. 

In this paper we have investigated electronic struc- 
ture and magnetic properties of Fei-a^Go^^Si alloys us- 
ing statical (LSDA) and dynamical (DMFT) mean-field 
approaches. The virtual crystal approximation gives 
us opportunity to investigate the electronic structure of 
Fei-a^Goa^Si system in whole range of concentrations. We 
have found that LSDA results strongly overestimate mag- 
netic moment values and extend magnetic phase dia- 
gram to much large values of Go concentration x com- 
paring with experiment. An account of correlation effects 
within DMFT results in good agreement with experiment 
for magnetic moment values as well as the position of 
the magnetization M{x) maximum. We came to con- 



clusion that Fei-a^Goa^Si is an itinerant electrons system 
which magnetic properties can be correctly described by 
LDA+DMFT method. 



A. Statical mean- field results 

The electronic structure of Fei-a^Go^^Si was cal- 
culated using the Tight Binding Linear-Mufiin-Tin- 
Orbital Atomic Sphere Approximation (TB-LMTO- 
ASA) method with conventional local-density approxi- 
mation (LDA).^^ The Fe and Go atoms in Fei-a^Go^^Si 
alloy were treated by virtual atoms with the atomic num- 
ber value averaged by the concentration x. The exper- 
imentally observed lattice constants^^ were used in the 
present band structure calculations. The calculated mag- 
netic moment as a function of cobalt concentration is 
presented in Fig. [11] and densities of states for x=0.2, 
0.5 and 0.8 are presented in Fig. [121 One can see that 
the calculated results start strongly deviate from experi- 
mental values of magnetic moments for Go concentration 
X > 0.3 with magnetization M{x) maximum at 0.5 jub 
instead of experimental value of about 0.2 /i^ and giving 
stable magnetism for larger x values than it is observed 
experimentally. Our results agree well with those pre- 




Concentration (x) 

FIG. 11: (Color online) Concentration dependences of the 
magnetization (in /is) obtained from neutron scattering 
experiment^ (blue dashed line), LDA+DMFT (orange bold 
line) and virtual crystal approximation (black thin line). 

sented in Ref|3l|l33. 

A small magnetic moment value calls the question of 
the nature of magnetism in Fei-^cGoa^Si: do we see itiner- 
ant electrons magnetism of the Stoner type or moments 
are local but their average value is suppressed due to 
quantum fluctuations and disorder effects? In order to 
investigate the localization degree of the magnetic mo- 
ment in Fei_a^Goa^Si system we have performed the su- 
percell calculation to simulate a Go impurity in FeSi. The 
supercell was constructed in the FeSi lattice with all ba- 
sis lattice vectors doubled and containing a total of 64 
atoms of transition metal. The impurity Go atom was 
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FIG. 12: (Color online) Density of states for x=0.2, 0.5 and 
0.8 obtained in LSD A calculations based on virtual crystal 
approximation. 



assumed to substitute one of Fe atom. The obtained 
magnetization spatial distribution in Fe63CoSi64 is pre- 
sented in Fig. [131 One can see that magnetic moment 
is not localized on the Co impurity. There is a mag- 
netic cluster containing one cobalt atom and six nearest 
iron atoms. This is a result of a strong hybridization 
of 3(i-states cobalt and iron through 35- and 3p-states 
of silicon. Such a magnetic cluster can be considered as 
a model for magnetic alloy Feo.84Coo.i6 Si for which the 
experimentally observed magnetic moment is about 0.13 
jiB per transition metal atom. The averaged value of the 
magnetic moment in the cluster of 0.13 /i^ is in excellent 
agreement with the experimental value. 

The fact that the magnetic moment induced in FeSi 
by Co alloying is delocalized supports the itinerant mag- 




0.003 



FIG. 13: (Color online) Schematic representation of a mag- 
netic cluster simulated with LSDA supercell calculations. 
Numbers are values of magnetic moments of atoms from dif- 
ferent coordination spheres. The arrows denote directions of 
the magnetic moments in the ground state. 



netism picture. In this case the Stoner approach is ex- 
pected to give a realistic description of Fei-a^Coa^Si mag- 
netic properties. However, the results of this section and 
previous works show that direct Stoner theory applica- 
tion leads to the overestimation of the magnetic moment 
value at intermediate concentrations (FigfTT]). This is a 
result of ignoring dynamical correlation effects that, as 
we will show below, result in a strong renormalization of 
states near the Fermi level and subsequent reduction of 
the magnetic moment value. 



B. Dynamical mean-field results 

In this section we investigate the influence of dy- 
namical correlation effects on the magnetic properties 
of Fei_a^Coa^Si. In our calculations we assume that 
Coulomb correlations treated by DMFT renormalize a 
paramagnetic density of states and then this new DOS is 
used as an input for Stoner theory calculations. 

Let us first qualitatively discuss the origin of mag- 
netism of the investigated alloys using the Stoner cri- 
terion 



IdNiEp) > 1. 



(8) 



Assuming the Stoner parameter 1^ value is equal to 1 eV 
the magnetic ground state is stable if the density of states 
at the Fermi level N{Ef) is larger than 1 states/eV. With 
Co substituting Fe in FeSi number of electrons per for- 
mula unit increases and the Fermi level runs through the 
peak above the energy gap. Hence for Fei-a^Co^^Si the 
value of N{Ef) is determined by the height of this peak. 
In DMFT the spectral function depends on the temper- 
ature. Fig. [T4l gives DMFT densities of states obtained 
by using the exact diagonalization technique. One can 
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see that the density of states value at the peak maxi- 
mum is larger than 1 states/eV at T= 58 K but becomes 
smaller than 1 at T=232 K. Then the Stoner criterion 
for magnetism Eq.® is satisfied for T= 58 K but not for 
T=232 K. This result agrees with experimental values of 
the Curie temperature (Tc = 50 K). 
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FIG. 15: (Color online) Density of states of the doped FeSi 
model obtained by using the exact diagonalization DMFT 
method at T=58 K for x=0.36 (brown solid line) and x= 0.66 
(red dashed line). The density of states near the Fermi level 
is presented in the inset. The black dashed line corresponds 
to the Fermi level. 



FIG. 14: (Color online) Density of states for FeSi obtained by 
using exact diagonalization DMFT method at T=232 K (red 
dashed line) and T= 58 K (brown solid line). 

We are now in a position to perform the quantitative 
analysis of the magnetism. The DMFT calculations for 
Fei_a;Coa;Si were performed for various concentrations x 
values. The results for x=0.36 and x=0.66 are presented 
in Fig. [151 One can see that for x=0.36 the density 
of states at the Fermi level N{Ef) > 1 and for x=0.66 
N{Ef) < 1. According to the Stoner criterion Eq.(|8]) 
that gives magnetic and nonmagnetic ground states for 
x=0.36 and x=0.66 correspondingly in good agreement 
with experimental data (see Fig. [TT]) . 

We have used the obtained paramagnetic DMFT densi- 
ties of states for different concentrations fFigfT5|) to solve 
the Stoner model. Self-consistent values for spin-up 
and spin-down numbers of electrons are given by equa- 
tions for the total magnetic moment 

/ + 00 
(iV(e + /dnt) - N{e + hn{))f{e, T)de, (9) 
-CO 

and the total number of electrons 

/ + CX) 
(N{e + hn^) + N{e + Idny))f{e, T)de (10) 
-OO 

that are recalculated iteratively. Here N{e) is a density 
of states obtained in DMFT calculations and f{e^T) is 
the Fermi distribution function. The Stoner parameter 
Id was chosen to be 1 eV close to that used in previous 
theoretical estimations.^^ 

The calculated concentration dependence of magneti- 
zation M{x) is presented in Fig. [Til There is good agree- 
ment between experimental and theoretical values. The 
main effect of using p{e) obtained in DMFT is a strong 
reduction of the resulting magnetic moment values com- 
pared to calculations using unrenormalized LDA DOS. 



This reduction is due to the quasiparticle weight factor 
Z = 1/ m* appearing in the numerator for Green function 
expression Eq.([5|). Then the integral over quasiparticle 
band states near the Fermi level is decreased by a factor 
of Z comparing with unrenormalized LDA values. As 
Z ^ 0.5 in our DMFT calculations that results in corre- 
sponding decrease of M{x) values by this factor. 



IV. DISCUSSION 

In this paper we have investigated electronic structure 
and magnetic properties of FeSi and Fei-a^Co^^Si systems 
using statical and dynamical mean- field approaches. Our 
band structure analysis supports the correlated band in- 
sulator model for this materials in contrast to the Kondo 
insulator model. The results of DMFT calculations have 
shown a strong renormalization of states near the Fermi 
level. The estimated band-mass renormalization m* ~ 2 
agrees with that obtained in the recent ARPES exper- 
iments. Analyzing paramagnetic DMFT densities of 
states calculated at different temperatures and at dif- 
ferent Co concentrations we have shown that itinerant 
magnetism picture is valid for Fei_a;Coa;Si alloys. 
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